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The renormalization-group properties of gauge-invariant transverse-momentum dependent 
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formalism and are absent when one uses a direct gauge contour. 
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I. INTRODUCTION 



Theoretical interest in the use of Wilson lines (also termed gauge links or eikonal phases) 
has been greatly stimulated recently by both theoretical investigations and experiments on 
single spin asymmetries (SSA) [l|, 0, 0, [s], 0] — see also for a quite recent review and 
further references. This interest appears in the context of a gauge-invariant formulation of 
parton distribution functions (PDF) in terms of hadronic matrix elements. Because these 
matrix elements involve quark-ant iquark field operators at different spacetime points, one 
has to introduce a path-ordered gauge link of the form 



[y,x\C] = Pexp 



y 

-ig I dz^A'^iz^ 

lx[C] 



that restores gauge invariance, albeit introducing an implicit, i.e., functional dependence on 
the contour C adopted. Here, C is, in general, an arbitrary path in Minkowski space and 
V denotes the path-ordering instruction that orders the Lie-algebra valued gluon fields with 
the earliest contour point x furthest to the right, whereas the final point y is put furthest 

(V) 

to the left. [Throughout this work, ta stands for the generator Ta of the fundamental 
representation (labelled F) of color SU(Ai'c). Note that there is an implied sum over the color 
index a.] 

The concept of gauge links is so pervasive throughout Yang-Mills theories because it 
ensures local gauge invariance independent of the particular dynamical theory. The concept 
of using contour-dependent operators in QCD is an old one and mostly studied in connection 
with the renormalization of singularities caused by contour obstructions, like end- or cross- 



points, and cusps 



21, 22, 23, 24, 25 



-see, for instance 



27 



@, i, m Q H El, 0, M, m E E, H M 

30, |3l| and further references cited therein. The renewed 



interest in Wilson lines in SSA and day-to-day applications is due to the potential breakdown 
of universality in transverse-momentum dependent (TMD) parton distribution functions (or 
PDF for short) js^, [H, H, S l36l . |37j and the ensuing discussion about a process dependence 
of the gauge link caused by the color flow of the particular hard-scattering process. As a 
consequence, it is argued that there is a change of the overall sign of the single transverse 
spin asymmetry from the deeply inelastic scattering (DIS) to the Drell-Yan (DY) case and, 
hence, a breakdown of universality j32|. The content and formulation of this question is 
coterminous with the properties of Wilson lines within the purview of specific gauge-fixing 
prescriptions and the imposition of boundary conditions at light-cone infinity. 

This brings in questions about the role of final (initial) state interactions (FSI) between 
the struck quark and the target spectators that may yield important effects like shadowing 
and SSA SHE. The core issue here is the choice of the gauge adopted for the Wilson 
line that has to be inserted to ensure gauge invariance. For singular gauges, like the light- 
cone gauge A'^ = 0, the gauge link vanishes by choice without any restriction put on the 
gauge potential at infinity. Then, to avoid light-cone singularities at = 0, specific 
boundary conditions on the gluon propagator have to be imposed in order to exhaust the 
remaining gauge freedom and recover the results obtained in non-singular gauges, say, in 

and further worked out and detailed in 



the Feynman gauge. It was shown in [41 



33 , that 



the effects of FSI in the light-cone gauge can be properly taken into account by including a 
gauge link (an "eikonalized" quark line) which involves a path in the transverse direction. In 
covariant, i.e., non-singular, gauges this additional term does not contribute and, therefore, 
the modified definition, proposed by Belitsky, Ji, and Yuan, for the TMD PDF reduces to 
the correct gauge-invariant one. These important findings not withstanding, yet a consistent 
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gauge-invariant picture for TMD PDFs in the whole phase space is still incomplete, because 
the behavior of the gauge contour at infinite transverse distance is largely arbitrary [4l|. The 
crucial point is — as we have recently shown [42] — that splitting the gauge link and allowing 
the separated contours to stretch out to infinity in the transverse configuration space, induces 
an additional contribution that cannot be dispensed with by imposing suitable boundary 
conditions [33|. Instead, one has to compensate this new contribution by incorporating into 
the definition of the TMD PDF an eikonal factor which provides a soft counter term in the 
sense of Collins and Hautmann 43|, 14^, |45 . 

In the present work, we shall investigate these issues from the point of view of the renor- 
malization group and address parton distribution functions — integrated and unintegrated — 
aiming for a more suitable definition of TMD PDFs. Our considerations will employ contour- 
dependent operators and we will calculate the leading gluon radiative corrections in the 
light-cone gauge. The renormalization of such operators is supremely simple to deal with 
when stated in terms of anomalous dimensions because these quantities are local and do 
not depend on the length of the gauge contour. Therefore, they provide a powerful tool to 
access the renormalization properties of Wilson lines and take into account those contribu- 
tions originating from geometrical obstructions, notably, endpoints, or sharp bends in the 
contour, as first pointed out by Polyakov jl]. In fact, we will show in more detailed form 
than in our brief presentation in (see also [1^1) that, adopting the light-cone gauge, the 
leading gluon radiative corrections associated with the transverse gauge link give rise to an 
anomalous dimension that exhibits a Inp"*" behavior — characteristic of a contour with a cusp. 
We will fathom out the physics underlying this finding — in particular, the renormalization 
effect on the junction point of two individual gauge contours joined non-smoothly through 
a cusp. Moreover, we will present a new definition for the TMD PDF that (i) reduces to 
the correct integrated case and (ii) coincides with the result obtained in the Feynman gauge 
that is untainted by contour obstructions (the reason being, we reiterate, that in this case 
A_L vanishes at infinity). 

The remainder of the paper is organized as follows. In the next section, we give a sum- 
mary of the kinematics used and sketch the spacetime picture of DIS. Section UTTl discusses 
integrated PDFs and their gauge-invariance and renormalization-group properties. In Sec. 
IIV[ we reinvigorate our statements on the transverse gauge link, presented briefly in [42], by 
a formal derivation, making use of a "classical" current along the lines of thought described 



by Jackiw, Kabat, and Ortiz in [47[. The calculation of the leading one-loop anomalous 
dimension of the TMD PDF in the light-cone gauge is outlined in Sec. |Vl Here we also 
present a generalized factorization rule for gauge links which takes into account the possibil- 
ity that the contours may be joined non-smoothly via a cusp-like junction point. The same 
section contains the evaluation of the soft counter term to supplement the definition of the 
TMD PDF and its interpretation as an "intrinsic Coulomb phase" , in analogy to the phase 
found by Jakob and Stefanis 48]for QED within a manifestly gauge-invariant formulation 



in terms of Mandelstam fields 49, M. In Sec. EH we turn our attention to the real-gluon 
contributions and the evolution equations, providing a tangible proof that the integrated 
PDF obtained from our modified TMD PDF definition coincides with the correct one with 
no artefact of the cusped contour left over. Section IVIII addresses the application of our 
approach to the Drell-Yan case. Finally, a summary and further discussion of our findings 
together with our conclusions is given in Sec. IVIIIi 
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II. KINEMATICS AND SPACETIME PICTURE OF (SI)DIS 

In what follows, we employ null-plane coordinates with 

= (P+^ P^) ^ p± = (po ^ p3^) 1^ ^ p2 ^ 2P+P- - Pi , (2) 



where P+ > is large and the other components P^, P_l > are small. Moreover, we will 
visualize the spacetime structure of the chief hadronic reactions, like DIS, Semi- Inclusive 
Deep-Inelastic lepton-nucleon Scattering (SIDIS), etc., as a series of snapshots on a plane of 
equal as depicted in Fig. [H [For a pedagogical exposition of this graphical method, see, 
e.g., Refs. 51, 13], and the review in Ref. |53].] 

Let us now fix the kinematics relevant for the cases to be considered in our work. We 
introduce two light-cone vectors 

n*'^ = l](l,l,Ox), n'^ = ^(l,-l,Ox) (3) 
with the following properties of their plus/minus light-cone components 

n*+ = V2Q , n*- = 0, n+ = , n~ = — ^ , n*n = 1 , {n*f = = , (4) 

where Q is an arbitrary parameter having the dimension of mass. Then, the momentum of 
the initial hadron reads 

pA* = + ^ p2 _ ^2 _ 

The momentum of the struck quark before being "measured" [lH by the photon is 

fcf^ = xP^" + k^^ k+ = xP+ = V2xQ (6) 

with 

Kn = , = (0+, k^) , fcin± = fcx , (7) 

whereas the off-shell photon has the momentum 



After the interaction with the highly virtual photon the struck quark acquires the momentum 

tt' M"^ + O"^ 

Ct = (C + = - ^>*' + 2x ■ 

Choosing a specific Lorentz frame upon setting Q = Q/{2x'), corresponds to an almost 
lightlike hadron that moves along the plus direction. In the Bjorken limit oo, the 

variables x' and x coincide up to 0{M'^ / Q'^) terms, so that one gets for the photon 

meaning that the photon moves along the negative x^-direction, whereas the struck quark 
(after being probed by the photon) moves along the minus direction to infinity: 




FIG. 1: (a) Schematic Feynman graph for DIS and (b) its spacetime picture at the amphtude 
level. The thick line denotes the struck quark, whereas the other solid lines along the direction 
represent spectators, and curly lines mark exchanged gluons. 



Therefore, in the Bjorken limit, one can estimate the corresponding conjugated spacetime 
light-cone coordinates to be 

5xp ~ 1/P+ ^ small , 5xj, ~ ^ large , (12) 

6x'i^ ~ — > large, ~ l/k" small . (13) 

In the next section, the spacetime picture, described above, will be used to motivate the 
introduction of the extra transverse gauge link. 



III. INTEGRATED PDFs: GAUGE-INVARIANCE AND RENORMALIZATION- 
GROUP PROPERTIES 



A well-known example for an integrated PDF is provided by the single parton distribution 
of a quark of fractional longitudinal momentum x and flavor i in a fast moving hadron which 
contains the nonperturbative physics in DIS:^ 

1 

2 1 27r 



fi/H{x) 



{H{P)\tPi{r,0^)j+UO-,0^)\H{P)) 



(14) 



where the quark momentum is defined by fc"*" = xP~^ with x = Q"^/ {2P ■ q). 

To give this expression a physical meaning, i.e., elevate it to an observable, one has to 
ensure that it is gauge invariant.^ To achieve this goal, one usually introduces ad hoc a 
gauge-link operator 



[y- 



X 



V exp 



r dz-A+ {O,z-,0^)ta 

J x~ 



(15) 



^ Strictly speaking, this is an unrenormalized quantity. The renormalization properties of the operators will 

be considered shortly in detail. 
^ Otherwise, one would have to define a "standard" gauge to be used in all calculations in order to make 

them comparable because what is a quark in one gauge is a quark plus gluons in another. 
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which is path ordered along a hghthke hne from the point x", where the quark was removed 
from the hadronic state by the annihilation operator ijji{0~ ,0±), to the point y' , where it 
was recreated by_0j(^~, 0_l). The gauge-hnk operator represents the struck quark as an 
eikonal parton |54|] with fixed color charge g, after its interaction with the highly virtual 

photon (0 < g2 = _g2 ^ I Q^y2^ ^j^g J3jg 

process, as it moves with (almost) the speed 
of light along a lightlike line in the x~ direction. The spacetime visualization of this process 
is illustrated in Fig. [I](b). 

It should be clear that, ultimately, expressions f|T^ and f|T5l) have to be quantized, using, 
for instance, functional-derivative techniques. This means that in employing Eq. f|T5l) in Eq. 
f fT^ . the gluon potential in the gauge link has to be Wick contracted with corresponding 
terms in the interaction Lagrangian, accompanying the Heisenberg fermion (quark) field 
operators. The consequence of this operation is that all terms in the expectation value in f|T^ . 
proportional to the gauge parameter and originating from the gluon propagator in a covariant 
gauge, will ultimately cancel. This has been explicitly proved for the quark propagator in 



leading-order perturbation theory long ago in [14'] (see also |l2l . Il3| and [55| for alternative 



formulations). Indeed, it has been shown in these works that after renormalization in a MS 
scheme, the anomalous dimension associated with the ordinary quark self-energy — which is 
gauge dependent — gets additional contributions, stemming from the two endpoints of the 



gauge link (coined in IJ] the "connector"),^ that exactly cancel its gauge-parameter term, 
so that the anomalous dimension of the composite gauge- independent quark propagator 
termed "hybrid" in [3]) is indeed free of the gauge parameter. [Below, the terminology of 



14l | is adopted.] This translates into the following sum rule for the anomalous dimensions 

Thybrid T2(7 ~l~ Tconnector ; (-^6) 

where 

7.. = ^ + 0{a!) , (17) 

7connector = '^^^f^^^ + 0{al) , (18) 

and 

7hybnd = + 0(«.^) , (19) 

with tts = 5'^/4vr, Cf = {N"^ — l){2Nc) = 4/3, and a being the gauge parameter. Note that 
we use the same conventions for the definition of the anomalous dimension as in ij], i.e., 
we write 

/i 1 (iZhybrid 



7hyb.d = ^^rr^ (20) 

^ ^hybrid "/^ 

with analogous expressions for the other anomalous dimensions. 

The above sum rule ( |T6l) is nothing but a "logarithmic", i.e., additive, version of the 
Ward- Takahashi/Slavnov- Taylor, identities in terms of ratios of the various renormalization 
constants of the QCD Lagrangian in the MS scheme 



■^connector-^lo Z3 Z3 NcO^s,^ , n1 /^i \ 

y ^ = ^ = ^ = 1 + (3 + a)-, (21) 

^hybrid ^1 Zi oTT 6 



^ These contributions are generated by singularities at the endpoints of the Une integrals that are dimen- 
sionally regularized in D = 4 — e dimensions and give 1/e poles. 
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where, ^hybrid = -^connector ■Z'2g, with ^25, ^3, Z3 being, respectively, the renormahzation con- 
stants of the quark, the gluon, and the ghost field, whereas Zi^, Zi, Z\ are the renormahza- 
tion constants pertaining to the quark-gluon-quark vertex, the three-gluon vertex, and the 
ghost-gluon-ghost vertex, respectively.^ 

Some important comments are here in order: (i) This leading-order result can be formally 
proved for smooth contours in every order of QCD perturbation theory (ii) It 

was shown in j^, [H, 13, 56] that the connector can be renormalized by multiplying it with an 
appropriate renormahzation constant and by replacing in the exponent the strong coupling 
by its renormalized version. Indeed, the renormahzation constants and Z^^ do not 
depend on the path chosen in any order of g (lol . 12 . [isl . 56] and the crucial Slavnov- Taylor 



identities are satisfied. This has been explicitly shown in [9|, |l3j up to the order g'^ and to 
all orders in g in [lo|. (iii) The straight line is actually enough for the renormahzation of 
the connector for any smooth contour p 



57 . Ssl . The reason is that what matters are 



only the singularities induced by the endpoints of the contour, which are multiplicatively 
renormalizable using exclusively dimensional regularization [13], whereas the specific path 
itself (for instance, its length) is irrelevant — provided no local obstructions like cusps and self- 
intersections are involved. Such obstructions would induce additional anomalous dimensions 
because of discontinuities in the contour slope, as discussed in detail in [sl, [H, 1^, that have 
to be taken into account in the anomalous- dimensions sum rule.^ We will show below the 
key role of a cusped gauge contour in the eikonalized TMD PDF with a transverse gauge 
link. 

The "eikonalized" quark PDF reads [s^] (see also 53, 60]) 



271 



.e-^'=^«-(P]^,(r,Ox)7+[r,0-]V^.(0-,Ox)]P) 



(22) 



and is a manifestly gauge invariant quantity, but has an anomalous dimension that comprises 
contributions stemming from the two endpoints of the lightlike contour Q-,o- (recall the 
remarks on the anomalous dimensions given above). Alternatively, one may be tempted to 
split the connector [^~,0~] into two gauge links that connect the points 0~ and ^~ through 
a point at infinity, the aim being to associate each of them with a quark field operator. 
Inserting a complete set of intermediate states, one can then recast Eq. (122]) in the following 
form 



split , 



rspi 
J i/a 



X] 



271 

~2^ 



'^^«-(P]^,(r,Ox)[r,oo-]t ]n)7+(n][oo-,0-]^,(0-,0x)]P) 



e-^+«"(p|fr^ (r,0±]Ci) ]n)7+(n]^, (0',Ox]C2) ]P) 



(23) 



where we have introduced the path-dependent Mandelstam fields 49, 50l 



^(xICo 



Pexp 



oo[C2 



(24) 



It is worth noting that in the covariant gauge a = —3, the gauge-contour divergences cancel among 
themselves and the Slavnnv- Taylor identity pTjl becomes trivially satisfied 
^ Just recently, Pobylitsa [59| has studied inequalities of a particular class of anomalous dimensions depend- 
ing on cusp angles. 
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FIG. 2: Space-time picture of the DIS process with a gauge connector (double hue) (a) and a spht- 
gauge hnk (two double lines) at light-cone infinity (b). The infinitely distant parts of the contours, 
where the behavior of the fields is not specified, are represented by the symbol for "ground" (earth) 
in electricity, introduced in 43] • Else, the same designations as in Fig. [T]are used. 



for the fermion and 



\l/(x|Ci) = 'ip^(x)Vexi) 



oo[Cx 
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(25) 



for the antifermion. These field operators represent the struck quark as an eikonal line 53l 
along a light-like ray in the minus light-cone direction while interacting with the gluon field 



51 



of the hadron, thus mimicking the motion of a struck quark in a real experiment 

The above two definitions ((Hj) and (l23l) are equivalent, because the anomalous dimension 
of the direct gauge link [^^,0~] is preserved when splitting the contour into two distinct 
contours through infinity. This means, in particular, that the junction point — which is shifted 
to infinity — is not creating any anomalous-dimension artefact. By virtue of the smooth 
connection of the contours Ci and C2, direct calculation shows that the renormalization of 
the junction point z preserves the algebraic identity 

[x2, z I Ci] [z, xi I C2] = [x2, xi I C = Ci U C2] (26) 

in any order of the coupling [l^ . This is not trivial because in a local gauge-invariant theory 
the factorized gauge link does depend, in general, on the junction point. However, for purely 
lightlike smooth contours Ci and C2, the above gauge-link factorization property is satisfied 
and the transport of color information by two different routes does not depend on the junction 
point (the latter being "hidden" at infinity in our case), with complete cancellation of the 
contributions at infinity. This is a well-established property which has been demonstrated by 
many authors in the early literature on studies of path-ordered exponentials. The analogous 
situation for non-smooth contours, which are not purely lightlike, and the generalization of 
Eq. fl26|) will be discussed in the next section. 



IV. DERIVATION OF THE TRANSVERSE GAUGE LINK 

To substantiate the use of the transverse gauge link and prepare the ground for the 
gauge-invariant formulation of unintegrated PDFs, which bear a full transverse-momentum 
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dependence, let us introduce a Coulomb source in terms of a "classical" current and write 

J,{y) = 9 fdy', 5^'\y - y') , y'^ = v^T , (27) 



which corresponds to a charged point-like particle moving with the four-velocity along 
the straight line f^t. The gauge field related to such a current has the form 



where D^^ is the gluon Green's function in an arbitrary covariant gauge. 
Here and in below, we will use two dimensionless light-cone vectors 



(2J 



n 



V2 



:i,o,,±ii 



[n 



[n 



-\2 



0, 



in^ ■ n 



and the metric tensor 

= + {n+Y{n-Y + {n+Y{n-Y with gi^" = -5^" 



(29) 



(30) 



Appealing to the spacetime structure of this process, illustrated in Fig. [2]J^a), we recast 
the current in the form 



n 



dr 6^^\y - n+r) + n~ / dr 5^^\y 
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5(2) ( 



n+5{y- 



dq e" 



iq 



271 q- + iO 



-n^5{y^ 



dq~^ e~ 



2tt g^ 



iO 



(31) 



which makes it clear that the first term in this expression corresponds to a gauge field created 
by a source moving from minus infinity to the origin in the plus light-cone direction, before 
being struck by the photon, whereas the second term corresponds to a gauge field being 
created by a source moving, after the collision, from the origin to plus infinity along a light- 
cone ray in the minus light-cone direction. Note that the underlying kinematics are those of 
Sec. nil 

To continue, we approximate the gluon Green's function D'^'^ by the free gluon propagator 
in the light-cone gauge = and obtain 



d^k 

Wy 



d^k 



-ik-: 



(27r)4 P + iO 



k^'in-Y + k^in^Y 



(32) 



where denotes the regularization prescription to handle the light-cone (pole) singu- 

larity at k^ . Taking into account that the n~-part of the current 0311) does not contribute in 
the = {A ■ n~) =0 gauge, one is able to obtain the transverse components /i = i = 1,2 
of the gauge field; viz.. 



-9 "K 



-^9 



^-^k■i ^/..(^) j dy+dy~d^y^e^^ y5{y~)5^^\yY] 



d^k 

2{2Tiy 

dk+ e-'''^^~ rd'^ki k\ 



27r [k+] J (27r)2fc2 



(33) 



The second, purely transverse, integral over d^fcj 
yields 



in Eq. (133|) gets factorized and finally 



2ti 



(34) 



(27r)2 kl 

while the first integral over dk~^ can be performed on account of the preferred regularization 
prescription for the pole at = 0. Note that in Eq. fl5^ . A is an auxiliary infrared (IR) 
regulator which ultimately drops out from all physical quantities. The connection between 
the regularization procedure in the momentum-space representation with the behavior of the 
gauge field at light-cone infinity can be anticipated from the following expression 



dk+ e-^*^+«" 



27r k+TiO 



IK 



dl 



dk+ 
( 

2-K 



dl 5{kI + r 



(35) 



where k = ±1. Then, one readily finds 



OO \K 



-V. 



An 



Coo Vln 



(36) 



emphasizing that this expression is dependent on the boundary conditions via the parameter 
Coo- Obviously, the longitudinal components vanish. 

Note that the magnitude of in our exp ression fl36l) turns out to be two times smaller 
than that obtained in Ref. [47|] (see also 6l| and references cited therein). The reason for 
this difference lies in the fact that in our case the source travels not along the whole plus 
lightlike axis, but only along half of it, changing its direction at the origin, as a result of its 
collision with the hard photon, and in agreement with the physical picture of the process. 

Let us continue by supplying within the same approach the gauge field in a general 
covariant gauge, characterized by the gauge parameter a, and using the gluon propagator 



b^'^ik) 



g^^ - (1 - a) 



Starting from Eq. (128|) . one gets after some simple algebraic manipulations 

A1 = 0, A'- = 0, A'+(0 = -^5(r)ln /|^^| , 

47r 



(37) 



(38) 



where fields in a covariant gauge are marked by a prime accent in order to distinguish them 
from those in the light-cone gauge [cf. Eq. (9) in Ref. 47|]. Notice that the a— dependent 
terms, which are proportional to ~ k~ under the d'^k- integral, do not contribute by virtue of 
the delta-function 5{k~). Next we give the (singular) gauge transformation which connects 
these two representations in the light-cone gauge and in covariant gauges: 



A 



LC 



0, m 



(39) 



Now we are ready to discuss the origin of the transverse contribution in Mandelstam's 
gauge-invariant formalism in the light-cone gauge. First, note that the analogous expression 
to the Mandelstam field fl2^ in a covariant gauge reads 



^cov(e;^+) =^exp 



-tg 



(40) 
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where the gauge field differs from the special case, A'~^, given by Eq. (138|) . Second, 
performing a regular gauge transformation 



exp 



dz'A^ 



(41) 



in the light-cone gauge on both sides of fl40|) . one can eliminate the Wilson-line integral in 
the phase. However, the regular transformation U{x~) does not exhaust the gauge freedom 
in the light-cone gauge completely and is, therefore, insufficient to trivialize the interaction 
of the struck quark with the gluon field of the spectators. More explicitly, a residual singular 
transformation f/sing(oo~, ^_|_) is still allowed and one realizes that the singular gauge trans- 
formation (!39l) refiects exactly this remaining gauge freedom. Carrying out this additional 
gauge transformation, one generates an additional phase that is now associated with the 
quark field itself; viz.. 



^(^ '^±)lC = f^sing(00 ,^x)V'Lc(^ 

l-ig rfz-4+,,,,(z-,|^) + 0(/) 



v^Lc(r,^±) , 



(42) 



which is now completely gauge-fixed, and hence, represents a quark with a fixed color charge. 
[This is indicated by a wide hat over the label LC which abbreviates 'light cone'.] Finally, by 
taking into account expression fl5Bl) . one finds that the quark wave function in the light-cone 
gauge acquires a phase that may formally be written as 



LC 



^Lc(r,i 



(43) 



The above arguments make it clear that a complete gauge fixing can be achieved in fl42l) 
by inserting the additional singular gauge transformation U^m^ which contains the cross-talk 
effects of the struck parton with the light-cone source. As a result, taking the product of 
two (local) quark field operators in the fixed light-cone gauge differs from what one finds in 
a covariant gauge. This difference is encapsulated in two phase factors so that one gets 



[^(r,^±)7+V'(o-,ox 



LC 



^Lc(^ ,|^)Pexp 



+ig 



dz I A. 



LC 



exp 



dz 



(ryO- 

^sourcel^ 



Ox 



(00 ,Z_ 
^LC(0- 



7^ 



(44) 



In the next section, we will use the results obtained above in order to demonstrate the role 
of the transverse link in the restoration of the prescription-independence of the anomalous 
dimension of the TMD PDF. Specifically, the explicit expression for the transverse gauge 
field at infinity, Eq. (1551) . will be used in the diagrammatic calculations of the gluon radiative 
corrections pertaining to the anomalous dimension of the TMD PDF. 



V. CALCULATION OF THE ONE-LOOP ANOMALOUS DIMENSION OF THE 
TMD PDF IN THE LIGHT-CONE GAUGE 

We have stressed before the importance of the renormalization effect on the junction point 
of the decomposed transverse contours at infinity. In this section we will explain exactly what 
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this means in mathematical detail. We will prove that the factorization of the gauge link into 
factors, each associated with a distinct contour starting (ending) at light-cone infinity, has 
to be modified to include an additional phase factor which accounts for the cusp anomalous 
dimension induced by the junction point of these decomposed contours. 



A. Definitions 



Our starting point is the operator definition of the (unpolarized) TMD distribution of a 
quark with momentum fc^ = {k~^, k~ , k±) in a quark with momentum = {p'^,p~ , 0_l): 



1 fd^-d^^s. 
2n{27i) 



(g(p)|V^(r,^x)r,^^;oo-,^Jt 



2 J 27r(27r)2 

X [c>o^,^_l; oo^, cx5j_]^7^[oo^, ooj_; oo^, 0_l] [oo^, 0_l; 0^, 0_l] 
where the lightlike and transverse gauge links are defined, respectively, by 



(45) 



[oo , z±; z , z±] = V exp 
[oo~, oo_l; oo~, = Pexp 



tg I dr A'^Xiz + n r) 





/■oo 

/ drl- AX{i±. + Ir) 
Jo 



(46) 



Let us emphasize that in the definition of the transverse gauge link the contour is defined in 
terms of the two-dimensional vector I which is absolutely arbitrary. We will show explicitly 
that this arbitrariness does not affect the local properties of the gauge link — in particular 
the anomalous dimension. 

Employing the light-cone axial gauge 



A-' 



-\2 



(A-n-) = , (n-) 







(47) 



the gluon propagator has an additional pole, related to the plus light-cone component of the 
gluon momentum, and reads 



D^C(n) = 



(48) 



To give this expression a mathematical meaning, we apply the following pole prescription 



1 



Ret/ Ad V 



g+ ± 17] 



PV 



1 1 

+ 



g+ + irj q 



+ _ 



(49) 



where t] has the dimension of mass, to be kept small but finite, and where we used the 
abbreviations Ret for retarded, Adv for advanced, and PV for the principal value. ^ In what 
follows, we regularize coUinear poles by means of the quark virtuality p"^ < 0, whereas IR 



We remark that another possible prescription — the so-caUed Mandelstam-Leibbrandt pole prescription 



62 



63 



6J| — is outside the scope of the present investigation, though we will make some related comments 



in connection with the anomalous dimension in Eq. (j79p . 
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+ {h.c.) 



FIG. 3: One- loop gluon contributions (curly lines) to the UV-divergences of the TMD PDF in 
a general covariant gauge. Double lines denote gauge links. Diagrams (b) and (c) are absent in 
the light-cone gauge. The omitted Hermitian conjugate diagrams are symbolically abbreviated by 
(h.c). 



singularities are regularized by an auxiliary gluon mass A which is put at the end back to 
zero. The described regularization procedure works well in the one-loop order, while at a 
higher loop order one may need to apply more sophisticated methods. 

In the tree approximation, where the gauge links are equal to unity and the quark-gluon 
interactions vanish, one trivially gets 

^ (50) 

= 6{p^ — xp'^)6'''^\k ±) -u{p)'y~^ u{p) , 

where u{p) denotes a quark spinor and summation over spin indices is tacitly assumed. 
Moreover, we use the short-hand notation 

u{p)-f+u{p) _ , . . 

^ = Mp) . (51) 

which implies 

f!il{x,^±) = 6(1 - x)S^'\k^) Mp) ■ (52) 



B. One-loop calculation 

Our objective here is to discuss the leading-order (LO) quark-gluon interactions, stem- 
ming, one one hand, from the standard QCD vertices, and, on the other hand, from the 
interactions of the quarks with the gauge links. Because in the TMD case the distance 
between the quark fields is spacelike, i.e., ^^^'^ = — ^ 0, the UV-divergent contributions 
arise only due to virtual gluon corrections. Therefore, in LO (alias, at the one-gluon ex- 
change level), diagrams (a), (b), and (c) contribute in a covariant gauge (see Fig. [3]), while 
in the light-cone gauge, only diagrams (a) and {d) give non-vanishing contributions — with 
the latter diagram being associated with the transverse gauge link. 

The Hermit ian-conjugate {h.c. for short) contributions (omitted in Fig. [3]) are generated 
by the corresponding "mirror" diagrams. In what follows, we consider first the "left" set of 
diagrams (in order to show explicitly how the transverse gauge link comes into play) and 
then we take the sum "left + right" which gives the total contribution to the TMD PDF. 
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These diagrams yield contributions proportional to delta-functions: 

where results from the diagrams shown in Fig. [31 The quark self-energy diagram (a) 

gives (in dimensional regularization with uj = 4 — 2e) 



S^^^P, as] ?7, e) = -g'^Cpfi 



2e 



with 



(27r)^ (p-g)2(g2_ A2 + iO) 
qf{n-y q''{n-Y 



p2 



(53) 



where the dependence on the auxiliary mass scale rj is "hidden" in the pole prescription [q 
and a = (7 ■ a) . 

The (7''''-proportional term gives a "Feynman" -like contribution, namely, 



(54) 



2e 



(27r)^ (p-g)2(g2 _ A2)p2 



(55) 



and generates no extra light-cone singularities. After carrying out the momentum integral, 
one gets 



Performing the remaining integral, one finally finds 

4eynman(P.as,/^,e) = -^CFr(e)(l - e) (^AtI 



xp'^ 



(56) 



p2 

l + e(2 + ^ln^!^-ln^^)+0(e^) 



(57) 

Note that the "mirror" diagram gives precisely the same contribution, doubling this result. 



C. Evaluation of the pole-prescription-dependent contributions 



The calculation of the [g"*"] -dependent part 

^oieb' ^; e) = g'^Cpfi^' 



»^ (p-g)2(g2 _ A2) 



— g)7+ 7^(]3 — q)q\ ip 



(58) 

is more demanding, owing to the presence of light-cone singularities, and we will consider 
its evaluation in detail. After some simple transformations of the numerator, we find 



''pole 



2e 



(Tq 
(2^ 



(rp g)^ 



1 



ip 



(g2- A2)[g+]p2 • 



(59) 



{p - g)2 

One observes that the integral / (^^2^x\q+] vanishes (which is true for the considered pole 
prescriptions but not for the Mandelstam-Leibbrandt one), while the rest can be recast in 
the form 



(60) 
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where 



dx 



'l-x) 



xil — x) \ 1 



(61) 



bearing in mind that by virtue of 7+7^ 
Thus, one has 



[n 



the term (72 (p) does not contribute. 



y(a) I y(a) 

Feynman pole 



,2\ « 



(P, ^, e) = ^Cp (^-47r^ ) r(6) <i (1 - 6) 



X 



X2_ 2 2_X2 

1 + e 2 + — In - In ^ 



A2 



2^"^ ^^.xii^n-eln 



+ 0(6^ 



X 



x(l — x) i 1 



xp'^ 



(62) 



In order to evaluate the integral J dx{l — x)/[x], one has to use a specific pole prescription 
for [x]. Let us consider three possible prescriptions: Advanced, Retarded and Principal 
Value: 



1 



1 



1 



Ret 



X + irj 



FJ Adv ^-^v 



1 



1 



PV 



1 



1 



2 \x + irj X — irj 



(63) 



using temporarily for convenience the short-hand notation fj = rj/p^ . In the limit of small 
r^, we keep only logarithmic terms and omit any powers of f/. The UV-divergent part (in the 
MS-scheme) then reads 



1 

47r e 



-— Of- 



1 A 27+p ■, V in . ^ 

1 - ln47r + 7s - ^-^ M + In ^ - — - m 

p+ \ p+ 2 



(64) 



where the numerical constant Coo depends on the pole prescription according to (cf. |33| ) 

C^= ( -1 , Retarded (65) 



, Advanced 

— 1 , Retarded 

— I , Principal Value . 



One the other hand, the finite part of the pole-prescription dependent gluon radiative cor- 
rections is 



''finite 



(p, a„ /i, ?7, e) 



a. 



27+p 



Cf 1 + In ^- + — In 



p2 p2 

7] in 



A2 



-In 



p^-X' 



V {.\ P 
dx- — ^ In 

m 



1 + In — ; — in Coo 1 In 



+ 2 
x(l — x) { 1 



p2 



xp^ 



(66) 



Evaluating this UV-finite integral ( 166|) by setting A^ = (which is justified given that 
this integral is IR finite and the magnitude of the IR regulator is irrelevant), we obtain 



1 (l-x) 

dx- In fx(l - x)l = 2 + (1 =F iri) 

xTiv 



Li, I - 

IT] 



Li2 



1 — if) 



(67) 
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As a result, the complete prescription-dependent finite part for = becomes 



Sfi"nite(P, V, e) = -TZ^F I 1 + In ^ 



a. 



fx- 



+2 + (iTiv) 



p. 
Li2 



p+ 



V 



1 + In — - inC^ In — 



Li2 



1 



1 =1= ir] 



(68) 



A key remark here is that any dependence on the pole prescription in Eqs. fl64p and 
is cancelled by taking into account analogous contributions originating from the transverse 
gauge link as depicted in Fig. [3](d). 



D. Contribution of the transverse gauge link at light-cone infinity 



The issue at stake in this subsection is the proof of the cancellation of all pole-prescription 
dependent terms in dealing with the gauge-invariant formulation of the TMD PDF. To be 
more specific, we focus on the effects related to the interactions with the gluon field of the 
transverse gauge link 



Pexp 



/"OO 

+ig (ir/^ ■ A_L(oo",0+;Z^r + |j 
Jo 



dTl±-A±{oo ,0~^;1±t) 



(69) 

and show that the contribution of the diagram [3](d) cancels out all terms proportional to 
the numerical factor Coo in Eqs. flMl) and fl66|) (or equivalently fl68|) ). Before we proceed, 
note that, as it is obvious from Eq. flMl) . the UV-divergent part depends not only on the 
pole prescription but also on the logarithmic p^-term. The effects related to this latter 
dependence will be considered subsequently. 

In Sec. IIV[ we have worked out the transverse components of the gluon field in the light- 
cone gauge and found Eq. (133|) . This expression can be further evaluated to read 



A^(oo-,0+; Z^r) 
finally assuming the form 

drl^ ■ A_l(oo", 0+; 1±t) = 



J 2tt J (27r)2 ^^^^ ' 



dq+ 



d'q± 
(2vr)2 



(70) 



(71) 



Consider now the free gluon propagator resulting from the correlation between longitudinal 
and transverse gluons; viz.. 



and use the relation 



q n 



(g2-A2)[g+] 



-0(27r)V^)(g + g') 



2mC^6{q+) 



(72) 



(73) 



to find the contribution of the diagram in Fig. [3](d): 



d'^q 

2\ « 



5{q 



7+(]3-g) 



iC^asC^ [ -47r^ ) r(e 



{p — qY{q'^ — P) 
(1 — x)6{x) 



2p+ Jq [x{1 - x)Y 



^1 ■ (^^' 



16 





FIG. 4: Graphical representation of a generic TMD PDF (shaded oval) in coordinate space. The 
double lines denote the lightlike and transverse gauge links, connecting the quark field points 
(0~,0^) and {£,~,^±), by a composite contour through light-cone infinity. The latter is marked 
by the typical symbol for the ground in an electrical circuit. The contour obstruction at infinite 
transverse and lightlike distance (oo~,oo_l) is symbolized by a cross, whereas the broken line 
indicates that this obstruction is "hidden" . 



One sees explicitly that by virtue of the relation 

1 



1 



lim 

v^o X ±ir] 



PV- T iTr6(x) 

X 



(75) 



the transverse-gauge link contribution f[74j) exactly cancels the dependence on the pole pre- 
scription in both the UV-divergent part S^y and in the finite part sjj^^j^g. To show this 
explicitly, we collect all pole-prescription dependent terms of diagram (a) in Fig. [3] and add 
to them the contribution from the transverse gauge link, i.e., Eq. (17^ . Then, we have 



lim / dx(l 



X] 



X — if] X + ir] 
-1 



lim / dx ln(x(l 

11^0 Jo X — ir] 



i2nCooS{x 
x)) 



In {x{l — x)) 



(76) 



which establishes the independence of the result on the parameter Coo — the latter encoding 
the adopted pole-prescription. As a result, the complete UV-divergent part of the TMD 
PDF fg/g{x,k^) is 



a.. 



{p,lj,,as;e) = Cp- 



TT 

Of - 

IT e 



Next, taking into account that 



1 
4 

1 - 



TP 
2p+ 

2p+ 



77 %7T 

1 + In — - - ivr Coo + ^vrCoo 



77 tn 

1 + n — 

2 



(77) 



2p 



+ 
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and recalling that the mirror (which we termed before "right" ) counterparts of the evaluated 
diagrams yield the complex-conjugated contributions, one can conclude that the imaginary 
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terms above mutually cancel, so that the UV- divergent part of diagrams (a) and (d) contains 
only contributions due to the p"*"- dependent term, notably. 



TT 



' ^ + ln^ ) -7i5 + ln47r 



(78) 



plus those terms originating from the standard MS renormalization. Hence, there is an extra 
anomalous dimension associated with the ^'''-dependent term which at the one-loop level, 
considered here, is given by 

7tSoop = [l + ^) = ~ ■ ^^^^ 

The difference S'y between 7smooth and 7lc is exactly that term induced by the additional 
divergence which ultimately has to be compensated by a suitable redefinition of the TMD 
PDF, if we want to reproduce the same anomalous dimension as in a covariant gauge. 
Here, some comments are in order. It was shown (see, e.g., (g^) that the Mandelstam- 



Leibbrand (ML) prescription [62|, |63|, |64 



' (80) 



[g+]ML g+ + ^0 q~ q+q- + zO 
yields a ^'''-independent anomalous dimension of the quark fields, i.e., 

l'^i. = l-C, + 0{al). (81) 

4 TT 

Moreover, the ML-prescription, in contrast to the PV/Adv/Ret ones (cf. Eq. 0651] ). entails 
additional poles in the complex g°-plane which allow for a Wick rotation and, therefore, 
it does not break the standard power counting rules. However, on the one hand, it is not 
clear how the ML-prescription can be related to any boundary conditions on the gauge field 
at light-cone infinity, thus making the popular initial/final state interactions interpretation 
questionable. On the other hand, the ML-regularization appears to be not sufficient for the 
calculation of the transverse gauge field at light-cone infinity (in the form of an expression 
analogous to, say, Eq. (H9l) ). The latter issue is potentially crucial for reproducing the results 
obtained in covariant gauges, while within the PV/Adv/Ret methods, the similarity between 
the light-cone and covariant gauges can be explicitly established. These issues will be further 
investigated and quantified in a separate work. 

To continue, recall that in a covariant gauge the gluon field vanishes at infinity and, hence, 
the only anomalous dimensions ensuing from the gauge link stem from its endpoints that 
are joined by a smooth direct contour.'' Actually, = [p ■ n~) coshx defines an angle x 
between the direction of the quark momentum p^ and the lightlike vector n~ . Then, in the 
large x limit, one has Inp"*" — > x, x ^ cxo- Thus, we come to recognize that the "defect" 
of the anomalous dimension, ^7, can be identified with the well-known cusp anomalous 
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7cusp(«s,x) = — Cf (xcothx - 1) , 

d , d a (82) 

— — - 57 = hm — 7cusp(as, X) = —Cf ■ 



It is worth reiterating that aU smooth gauge contours yield the same anomalous dimensions, 7smooth, as 
the straight line between the endpoints, because only the latter are relevant. 
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FIG. 5: Renormalization effect on the junction point due to gluon corrections (illustrated by a 
shaded oval with gluon lines attached to it) for (a) two smoothly joined gauge contours Ci and 
C2 at point 3 and (b) the same for two contours joined by a cusp (indicated by the symbol (8)) at 
infinite transverse distance (marked by the earth symbol) off the light cone. All contours shown 
are assumed to be arbitrary non-lightlike paths in Minkowski space. 

This is an important observation that deserves to be discussed in some detail. 

As we mentioned earlier, splitting the gauge contour for non purely lightlike contours 
through the light-cone infinity, is not equivalent to the situation with a direct contour between 
the two field points. To understand the deeper reason for this difference, we have to study 
again the algebraic identity (!26l) for decomposing (factorizing) gauge contours (links). The 
crucial question here is whether the defect of the anomalous dimension, we calculated, is 
compatible with this identity when the junction point is assumed to be at infinite distance 
in the transverse configuration space. To answer this question, consult Fig. [51 Panel (a) 
of this figure shows the renormalization effect (illustrated by a shaded oval with gluon lines 
attached to it) on the junction point in the algebraic identity (12^ . The contour Ci U C2 
is smooth and non-self-intersecting owing to the assumption that Ci and C2 are smoothly 
connected at 3. Then, both contours the direct one, C, and the decomposed one, Ci U C2, 
between the endpoints 1 and 2, cannot be distinguished from each other by switching on 
gluon quantum corrections. In particular, no anomalous dimension emerges from the junction 
point 3 and thus (symbolically) 

7c = 7C1UC2 ■ (83) 

Now we may ask what changes are induced, if we allow the junction point 3 to be shifted to 
infinity in the transverse direction off the light cone. The graphics at right of Fig. [5] helps 
the eye catch the key features of the situation involving two non-lightlike contours Ci and 
€2- It turns out that the naive assumption that 

7c = 7cf°uc|° (84) 

is incorrect for contours containing transverse segments. Instead, we found that in this case 

7c = 7cf=uc|° + 7cusp • (85) 

Consequently, the validity of the algebraic identity ( (|2B1) ) is not conserved and we have to 
replace it by the generalized gauge-link factorization rule 

[2,1\C] = [2, oo|Cr]^[oo, l|Cr]e^*-- , (86) 

which is valid for arbitrary paths in Minkowski space. In this expression, $cusp takes care 
of the effect induced by the cusp-like junction point. We will consider an explicit example 



19 



of such a phase in the next subsection, where we show how to compensate it by an eikonal 
factor in the definition of the TMD PDF. One may associate this phase with final (or initial) 
state interactions, as proposed by Ji and Yuan in |41i], and also by Belitsky, Ji, and Yuan 
in 3^. However, these authors (and also others) did not recognize that the junction point 
in the split contour (taking a detour to light-cone infinity) is no more a simple point, but 
becomes a cusp obstruction ~ Inp"*" that entails an anomalous dimension as the result of a 
non-trivial renormalization effect owing to gluon radiative corrections. 

These arguments make it clear that the naive decomposition of gauge contours that stretch 
out to light-cone infinity along the transverse direction is erroneous, simply because the basic 
algebraic identity fl26|) . which is tacitly assumed, is inapplicable to such contours and has to 
be replaced by Eq. fl86l) . It almost goes without saying that the modified factorization rule, 
expressed through Eqs. fl85|) and fl86l) . is valid when one is composing non- smoothly any gauge 
contours with a cusp obstruction at the junction point. What marks out a cusped contour 
from all the others, however, is that it gives rise to an anomalous dimension proportional to 
Inp"*", i.e., to a jump in the four-velocity. This is a salient ingredient in describing correctly 
a DIS process in spacetime, because if the two quarks (the struck one and a spectator) are 
separated also in the transverse coordinate space, the gluons emitted mismatch in rapidity 
and, hence, the contour liaising them has to have a sharp bend and cannot be the direct one. 



E. Compensating the defect of the anomalous dimension by a soft counter term 

The defect of the anomalous dimension owing to the gauge-contour cusp at light-cone in- 
finity represents a distortion of the gauge-invariant formulation of the TMD PDF in the light- 
cone gauge. To restore its consistency, we have to dispense with the anomalous-dimension 
artefact of the cusp. This can be achieved by supplying the original definition of fq/q{x, k±) 
by a soft counter term in the sense of Collins and Hautmann jH, 44, 45, 66| : 



i? = $(p+,n"|0)<l>"^(p+,n"|O , 
where the eikonal factors are given by 



^7) 







$(p+,n"|0) = 
$t(p+,n-|0 = (o 



V exp 

V exp 



ig / rfC^ f^AliC) 

Ccusp 



Ccusp 



(89) 



and evaluate R along the non-smooth (non-lightlike) integration contour Ccusp, defined by 
Ccusp : C = {b>, -oo < s < 0] U [n-s\ < s' < oo] U [l^r, < r < oo]} , (90) 

with n~ being the minus light-cone vector, as illustrated in Fig. [61 

Contour fl90l) is obviously cusped: at the origin, the four- velocity p+, which is parallel to 
the plus light-cone ray, is replaced — non-smoothly — by the four- velocity n~ , which is parallel 



A similar factorization rule holds for contours joined through a self-crossing point. 
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(0-. -O0+, Oi) 

FIG. 6: The integration contour associated with the additional soft counter term. 



to the minus light-cone ray. This jump in the four-velocity becomes visible in the standard 
leading-order term 



9 



ds I ds' 



{ViS - V2S')' 



9 



ds I ds' 



(f 1 ■ V2) 



vfs + f|s' — 2(t>i ■ V2)ss' 



(91) 



in which vi = p~^, V2 = n~ , and the change of the four- velocity at the origin produces an 
angle-dependence via (fi ■ ^2) = p"*". This means that exactly at this point the contour 
has a cusp that is characterized by the angle x ~ Inp"*" = ■ n~), and, therefore, the 
corresponding eikonal factor (!89|) gives rise to a cusp anomalous dimension. Obviously, this 
is exactly what we need in order to compensate the extra term in the anomalous dimension 
found in the preceding subsection. 

Next, we show that the one- loop gluon virtual corrections, contributing to the UV diver- 
gences of R and displayed in Fig. U\ yield an anomalous dimension that neutralizes the cusp 
artefact ^7. Note that in the light-cone gauge yl+ = (n^ ■ A) = only the first lightlike 
ray — 00 < s < and also the transverse segment contribute, since the other eikonal line 
along the minus lightlike ray depends on the longitudinal component of the gauge field and 
vanishes due to the gauge condition. 

Calculate first the diagram (a) in Fig. [71 In leading order, the first nontrivial term in Eq. 
(EID reads 



da / dr 



d^q e-"^- 



iq-u{a—T) 



q^'n"" + q^'n'^' 



d'^q 



(27r)^ g2 _ 



(27r)- g2__x2 
+ 



{q ■ u — z0)2 {q ■ u — iO) 



(92) 



The first term in the square bracket vanishes since is chosen to point along the p'^- 
direction, i.e., = (p^, 0^, 0^) , = 0, and by recalling that in dimensional regularization 
u^/{u^~'^) = 0. Notice that in a covariant gauge this diagram would be tantamount to 
the self-energy contribution of the struck quark. However, in the light-cone gauge, we are 
employing, a second term in the parenthesis in the first line of Eq. (p2l) appears which stems 
from the gluon propagator in that gauge. This term, being not lightlike, also entails a 
contribution to the cusp-dependent part, as we will now show. Indeed, one has 



<l.(^-^-P)(«,r/) 



d'^q 



{27r)^ (g2- A2)(g.M-iO)[g^ 



(93) 
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+ (h.c. 



FIG. 7: Virtual gluon contributions to the UV-divergences of the soft counter term, given by Eq. 
The designations are as in Fig. [3j 



an expression which would correspond to a vertex-like contribution of the pure gauge link in 
a covariant gauge — as one may appreciate from Eq. (!9T|) . 

The pole-prescription dependent integral can be evaluated in analogy to our previous 
calculations in the preceding subsection, so that 



[-a 



(94) 



where the bracketed term in the denominator is to be evaluated with the aid of Eq. (1751) . 
The last step in obtaining an explicit expression for ^"°^\u,r]) is to carry out the line 
integral over x (which enters because of the appearance of the gauge link in the lightlike 
direction). To do so, we have to take care of the additional (logarithmic) singularity ~ Inr, 
which cannot be regularized by the parameter 77, where r is an extra regulator. The origin of 
this singularity is related to the vector which defines the lightlike direction. There are, of 
course, several possibilities how to regularize this type of integral. For instance, one can get a 



regular expression in r — after the integration over drj — as discussed in Refs. ([19||, [2l|]) having 
recourse to the fact that the derivative (9$i^ ^"°^\u,ri)/dri is r-independent. However, for 
technical convenience, we apply here a different regularization technique by making use of an 
auxiliary regulator r and absorb the light-cone singularity In r inside a redefined parameter 
fj = 2t7], the latter being contained inside the pole-prescription contribution — cf. Eq. fITS]) . 
This is possible, given that r does not depend on the scale parameter rj and hence does not 
contribute to the evolution of the considered quantity (see Sec. rvT|) . 

Performing all these operations and taking into account that = P~^, one gets for the 
UV part of diagram (a) in Fig. [7] 

^V'^'^V) = --Cf- fin 4 - 4 - ^^^-1 • (95) 

Evidently, this expression yields a cusp-dependent contribution, ~ Inp"*", as we already 
mentioned, which would be completely absent in a covariant gauge, thus, underlying their 
mutual difference. Besides, there is a dependence on the choice of the pole prescription (via 
the numerical parameter Coo). In order to cancel this latter dependence, one needs to take 
into account the contribution of diagram (d) in Fig. [3 viz.. 



(27r)2(g^-A2)[g+] 
q ■ u + iO q'j_ ■ l± + iO 



x(-z)(27r)V^)(g + g') V, , ,^ • (96) 
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Using Eqs. (^-(^, we find 

$(d) = ig^CFfl^'TcC^ 



(27r)^-2 ql + P 



Li 

-asCYiTiCoo^{t) [ -47r— 



2\ « 



and extracting the UV-pole and adding it to Eq. (l95i) we finally arrive at 



(a+d) , 
UV 



TT 



[f]) = -Cp- ( In — — i «7rCoo + «vrC( 



TT 



pi 



TT 



pi 



(97) 



(98) 



This result exhibits the independence on the pole prescription of the soft factor R, in close 
analogy to Eq. (179ll . Taking into account the corresponding "mirror" diagram (which doubles 
the real part and cancels the imaginary one), we obtain the total UV-divergent part of the 
soft factor in the one-loop order: 



(1-loop), 
UV 



TT 



(99) 



making it apparent that there is no dependence on the pole prescription, as now all Coo~ 
dependent terms are absent. 

To conclude, we have shown at the one-loop level that the soft counter term (soft eikonal 
factor) has the following two important properties: 

(i) it gives rise to the same cusp anomalous dimension as fg/q{x, k±), but with an opposite 
sign, and 

(ii) it bears no dependence on the choice of the pole prescription to go around the light-cone 
singularity in the light-cone gauge (with corresponding terms cancelling among themselves). 

Therefore, it is reasonable to redefine the conventional TMD PDF by including into its 
definition the soft counter term ab initio. This provides 



(x, fcx; v) 



1 /'c/^_^^_^ 



X [oo^, oo^, oo_l]'''7+[oo^, oc±; oo^, Oj_] [c>o^, Oj_; 0^, 0_l] 
xij{0-, Oj_Mp)) h{p+, n-\0-, J$^(p+, n-\r, ^± 



(100) 

which represents one of the main results of our investigation here and in [42l |. 

Let us finish this section by gi ving a physical interpretation to the soft counter term, 
using Mandelstam's formalism [49|, l50[]. To this end, we utilize the exponentiation theorem 
for non-Abelian path-ordered exponentials [13] and recast the exponential operator (!89|l in 
the form 



$(m, n ) = exp 



n=l 



(101) 



where the functions have, in general, a complicated structure that is, however, irrelevant 
for our purposes here. The leading term in this series, $i, is just a non-Abelian generalization 
of the Abelian expression 



^ Ccusp 

Then, by virtue of the current 



(102) 



(103) 
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evaluated along the contour Ccusp (cf. Eq. (pOj) ) and where the velocity Vy equals either Uy, n~ , 
or Zj_ (depending on the segment of the contour along which the integration is performed), 
one can rewrite fll02p as follows 

^i{u,n~) = -fAn I dx^ j d^z5''''D'"'{x - z)jI{z) . (104) 

Ccusp ^ 

This expression looks formally very similar to the "intrinsic" Coulomb phase found by Jakob 
and Stefanis (JS) [11] in QED for Mandelstam charged fields involving a gauge contour 
which is a timelike straight line. The name "intrinsic" derives from the fact that this phase 
is different from zero even in the absence of external charge distributions. Its origin was 
ascribed by JS to the long-range interaction of the charged particle with its oppositely 
charged counterpart that was removed "behind the moon" after their primordial separation.^ 
This phase is acquired during the parallel transport of the charged field along a timelike 
straight line from infinity to the point of interaction with the photon field and is absent in 
the local approach, i.e., for local charged fields joined by a connector. It is different from zero 
only for Mandelstam fields with their own gauge contour attached to them and keeps track of 
its full history since its primordial creation. Keep in mind that the connector is introduced 
ad hoc in order to restore gauge invariance and is not part of the QCD Lagrangian. In 
contrast, when one associates a distinct contour with each quark field, one, actually, implies 
that these Mandelstam field variables should also enter the QCD Lagrangian (see [4g] for 
more details). However, a consistent formulation of such a theory for QCD is still lacking 
and not without complications of its own. 

The analogy to our case is the following. First, formally adopting a direct contour for the 
gauge-invariant formulation of the TMD PDF in the light-cone gauge (Figure [Ht^a) shows 
an example of the contributing diagrams), the connector gauge link does not contribute any 
anomalous dimension — except at the endpoints; this anomalous dimension being, however, 
irrelevant for the issue at stake. Hence, there is no intrinsic Coulomb phase in that case. 
Second, splitting the contour and associating each branch to a quark field, transforms it 
into a Mandelstam field and, as a result, adding together all gluon radiative corrections at 
the one-loop order, a ^'''-dependent term survives that gives rise to an additional anomalous 
dimension. We have shown that this extra anomalous dimension can be viewed as originating 
from a contour with a discontinuity in the four-velocity x{a) at light-cone infinity — a cusp 
obstruction. 

Classically, it is irrelevant how the two distinct contours Ci and C2 in Fig. [5] are joined, 
i.e., smoothly or by a sharp bend. But switching on gluon quantum corrections, the renor- 
malization effect on the junction point reveals that the contours are not smoothly connected, 
but go instead through a cusp. Here, we have a second analogy to the QED case discussed 
above. Similarly to the "particle behind the moon" , this cusp-like junction point is "hidden" 
and manifests itself only through the path-dependent phase fll04p . Note in the same context 
that integrating over the transverse momentum (see next section), the p^-dependent terms, 
resulting from virtual gluon corrections, cancel against their counterparts from real-gluon 
corrections, so that this cusp-induced phase disappears [see for illustration Fig. [8]J 



In our previous paper [42|, we concentrated on the anomalous dimension of the TMD PDF, 



and, therefore, only the UV-divergent parts were studied. In the present work, however. 



^ The existence of a balancing charge "behind the moon" was postulated before by several authors — see 
for related references — in an attempt to restore the Lorentz covariance of the charged sector of QED. 
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(a) (6) 

FIG. 8: (a) Schematic illustration of the direct gauge contour with virtual gluon-line insertions 
(denoted by curly lines), (b) Similar illustration for real gluon exchanges over the cut, the latter 
indicated by a dashed line. In both panels double lines represent the gauge links. 



we take also into account the UV-finite parts and, consequently, the dependence on the 
transverse momentum appears explicitly, as we discussed above. 



VI. REAL-GLUON CONTRIBUTIONS AND EVOLUTION EQUATIONS 

In this section, we concentrate our efforts on two subjects: (i) First we discuss in some 
detail the evolution behavior of the TMD PDF and establish the connection between our 
approach and that of Collins and Soper [g^]. (ii) Second, we prove that the integrated PDF, 
obtained from our modified definition, coincides with the standard one with no any artefact 
of the cusped contour used in the TMD PDF left over. 

(i) Evolution behavior. The modified TMD PDF f llOOl) depends on two arbitrary mass- 
scale parameters: the UV scale /i and the extra regulator r]. The yU-dependence is described 
by the standard renormalization-group evolution equation (see below) and is controlled by 
the UV-anomalous dimension, which arises as the sum of the anomalous dimensions of all 
the ingredients of the TMD PDF ffTOOl) : 



7/,/, = 72q + XI ^gauge link + 7R 
i=l 

j^C^ + 0{al). (105) 

4 71 



The anomalous dimensions associated with the quark fields and the soft counter term R are 
marked by self-explaining labels. We have used for convenience a short-hand notation to 
denote the anomalous dimension of each gauge link on the right-hand side of Eq. fllOOp by a 
number in the order the gauge link appears from the left to the right. 

Before we proceed, a couple of important remarks are here in order. One realizes that 
the anomalous dimension of fg/g coincides with the anomalous dimension of the conventional 
quark propagator in the light-cone gauge, but with the opposite sign due to the different 
Dirac structure. Up to the sign, this result also coincides with the anomalous dimension of 



the gauge-invariant quark propagator in a covariant gauge . The anomalous dimension of 
R, 7r, cancels precisely those contributions in the sum above which contain the ^'''-dependent 
terms. 
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FIG. 9: The leading-order real gluon contributions to the TMD PDF are shown. The diagrams (b) 
and (d) with a transverse gauge link do not contribute to the TMD PDF in the light-cone gauge. 
The dashed line marks the cut. 



On the other hand, the dependence on rj is more complicated and is described by an 
integral kernel to be determined below. In order to derive the corresponding evolution 
equation, one needs to calculate the real-gluon contributions, depicted in Fig. O Here, we 
present this calculation in the small-77 limit (that corresponds to the large-rapidity ( ^ 00 
limit within the CoUins-Soper approach [67]). Let us emphasize that in the case of the 
integrated PDFs, where the dependence on the regularization parameter rj appears at the 
intermediate steps of the calculations (in the light-cone gauge), it cancels out in the final 
expression. This will be demonstrated below. In contrast, in the unintegrated PDFs, this 
dependence remains and, thus, it should be treated by means of a corresponding evolution 
equation. Note that the diagrams in Fig. M are UV-finite and do not contribute to the 
anomalous dimensions. However, they do depend on the regularization parameter t]. The 
diagram (a) yields 

^real ' 9 J ^^^y _ DlSC [19 [q)\ 

x5{p-^ ~k+ -q+)5^^\q^-k^) , (106) 
where the absorptive part of the gluon propagator reads 

Disc [D^^^iq)] = 2ne{q^)6{q' - A^) i^-g^'' + j • (107) 

In the last equation we have adopted the PV-prescription, because the real gluon contri- 
butions are prescription-independent and the diagrams with transverse gauge links do not 



contribute. The ?7-divergences can be isolated by means of the standard rules given in [68 



^ ^ ^(1-a:) ln^+ ,^ \ . (108) 



After some standard calculations, one gets the "Feynman" (r^- independent) part 

y.(a)real _ ^ |1 - x\ k\ + x\'^ + x{x - 3)p^ 

^Feyn^an - [fc^^ + ^^A^ _ a:(l - x)j>f ^ ^ 

The 77-dependence appears through the pole-contributions, i.e.. 



y-i(a)real ^ 



^ - 5(1 -x) In 4 



(1 - x)+ 



+ xX^ — x{l — x)p'^ 



;iio) 
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On the other hand, the diagram (c) in Fig. [9] yields 

poo poo p M 

sSil = ^9' u,u, 6{p^ - xp+) da dr j ^^e-^'^M^-^)6('\q^ - k, 



qf^n + q n 
b+]pv 

In the small-?7 hmit, a straightforward calculation gives 



'1111 



"'real 



— Cf6{1-x)-^ 



2t\ 



1 - In- 



pi 



(112) 



Finally, the (logarithmic) dependence of the modified TMD PDF on ri is determined in terms 
of the equation 



TT 



Cf5(1-x) 



1 



1 



TT Vfcl +XA2 -X(l -X)p2 fcl + A2 



0o(p) • (113) 



We can recast this equation, which governs evolution with respect to in a form which 
formally resembles the standard CoUins-Soper evolution equation 67, 6§| with respect to /z, 
namely, 

v-^fj.'^ix, k^; yu, r]) = [/C(/i) + ^(/i, v)] ® f^fix, fcx; v) ■ (114) 

The renormalizat ion-group behavior of the functions /C(/i) and ^(/x, 77) [t^] is determined by 
the universal cusp anomalous dimension 



;ii5) 



Extracting explicit expressions for /C(/u) and Q{fi,ri) from our Eq. ( lllSp . we can readily 
show that they each satisfy Eq. (11151) with respect to the cusp anomalous dimension. We 
emphasize that the parameter 7] in our approach plays a role akin to the rapidity parameter ( 
in the additional evolution equation of CoUins and Soper, with Eq. (11131) being the analogue 
of the Collins-Soper equation. 

Therefore, the dependence on the dimensional regularization scale /i of the re-defined 
TMD PDF (llOOl) is given by the following renormalization-group equation 



(116) 



It is important to appreciate that only the modified TMD PDF, given by Eq. (llOOp . obeys 
such a simple UV-evolution; without the soft counter term, non-trivial extra contributions 
would arise in the corresponding anomalous dimension on the right-hand side of Eq. (llOSp . 
Taking logarithmic derivatives of /^°'^(x, /i, r^) with respect to both scales /i and 77, we 
get 



d 

^ d^ 







(117) 
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which estabhshes the formal analogy between our approach and the Colhns-Soper one. This 
equation ensures the absence of extra UV-singularities related to artefacts owing to the hght- 
cone gauge and is equivalent to our initial requirement of the cancellation of undesirable 
divergences. 

(ii) Integrated modified PDF. Consider now the integration of Eq. fllOOp over fc^. We 
collect all ?7-dependent terms from the virtual and the real-gluon contributions — contributing 
UV divergences — and perform the k± integration. We find that the result is r^- independent, 
so that the DGLAP evolution of this quantity is guaranteed. Below, we demonstrate this 
cancellation explicitly. The integration of the UV-divergent term fl78|) trivially gives 



2^CF-\nf]6{l-x) [d^k±6^^\k^) =2—CF-\n^6{l-x) . (118) 
vr e J TT e p+ 

On the other hand, the integration (in the dimensional regularization) of the T^-dependent 
real-gluon contribution ( IllOp yields 



^Cpln45(l-x)/i2^ d'-'^k^-^^^ =-2^Cpr(6) In 4 

p+ ' 'J ki+A^ TT p+ 



x6{l-x)(^^y , (119) 
where = xA^ — x{l — x)p'^. After extracting the UV-divergent term 

-2— Cf - \n^6(l-x) , (120) 

TT e p+ 

we observe that it exactly cancels the right- hand- side of Eq. flll8p . The same cancellation 
occurs between the r^-dependent terms in the virtual and the real-gluon contributions of the 
soft factor. Integration over the transverse momentum (using dimensional regularization) in 
the modified TMD PDF ffTUUD yields— at least formally— the integrated PDF 

J d--^kU^:\x, kr, /i, T]) = fi,a{x, /X) (121) 

which bears no 77-dependence as well. 

From the above considerations it becomes apparent that the renormalization-group prop- 
erties of this distribution are described by the DGLAP equation 

d ^ r (Ji^ / ^ \ 



where the integral kernel reads (in leading order) 



Pij{x) — 



-6{1 -x) + 



+ 0(«^) (123) 



2 ' ' [l-x) 
and the ()+-regularization is defined in the standard manner by 

Thus, one may conclude that the extra cusp-dependent terms are not present in the inte- 
grated case and, consequently, the UV properties of the standard PDFs (governed by the 
DGLAP equation) are not affected by the additional parameter rj, as expected. 
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(a) 



FIG. 10: (a) Spacetime picture of the Drell-Yan process at the amphtude level. The thick line 
denotes the struck quark, whereas the other solid lines along the direction represent spectators, 
and curly lines mark exchanged gluons attaching to the gauge links (double lines), (b) Splitting 
the gauge link in the DY process in analogy to SIDIS in Fig. [2] and using the same designations as 
there. Notice that the two gauge links are separated by a transverse distance off the light cone. 



VII. DRELL-YAN AND UNIVERSALITY 

The study presented in the preceding sections was performed for the semi-inclusive DIS 
(SIDIS). In this section we discuss the DY case within our approach and comment on uni- 
versality. 

We now ask ourselves to what extent our analysis can be applied to other reactions, 
like the DY lepton-pair production. It was shown by Collins (32j that the direction of the 
integration contours in the lightlike gauge links entering the definition of the TMD PDF 
should be reversed relative to the SIDIS, see Fig. [TOl 

[r,^±;oo",^j_]^[oo-,Or,r,o±Lu,s (125) 

[r,li;-oo-,^J[-oo-,0^;0-,Ojt^ . 

Later, Belitsky-Ji-Yuan [sl] argued that the transverse gauge links, they had introduced to 
exhaust the gauge invariance in the light-cone gauge, should be used in the DY case with 
the reverse sign, keeping in mind that the lightlike gauge links do not contribute. Hence, 
one has for DY the following combination of gauge links 

[oo~, oo^, oo_l]^'[cx)", oo_l; oo", Ojgjnis — ^ (126) 
[-00", -oo", oo_l] [-00", oo_l; -oo', OjJ^^ . 

In our approach, the latter replacement should be supplied with a change of sign of the 
additional regulator r] (cf. the pole prescription in Eq. (H9l)). i.e., 

^siDis — ' -VuY ■ (127) 

This change reflects, in fact, the different behavior of the gauge fields at the plus and minus 
light-cone infinity subject to the proviso of different pole-prescriptions. In the non-polarized 
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case (we exclusively discuss in the present investigation), this affects only the purely imagi- 
nary terms of the UV-divergent parts — consult Eq. Thus, in the DY case, one has 



= -— Gf- 

DY Arc e 



l-ln4vr + 7i.-^fl + ln^ + |^ + mC°Y 



pT y p-r 

where the numerical factor differs from the SIDIS case and is defined as 



:i28) 



DY 

oo 



-1 , Advanced 

, Retarded 

1 , Principal Value 



(129) 



These imaginary terms occur at the intermediate steps of the calculations, but do not con- 
tribute to the final expressions for the unpolarized TMD PDFs. Exactly the same arguments 
hold for the additional soft factor fl92|) . This means that the (real- valued) UV anomalous 
dimension of the unpolarized TMD PDFs is universal as regards the SIDIS and the DY 
processes: 

(130) 



SIDIS 

Ifq/q 



Ifq/q 



This, however, may not be true for the spin-dependent TMD PDFs, since in that case the 
imaginary parts play a crucial role and, thus, a sign change (expressed in (11281) ) might indeed 
affect the renormalization-group properties and the corresponding evolution equations. This 
is an interesting task which will be pursued separately elsewhere. 



VIII. SUMMARY AND CONCLUSIONS 

In this paper we have applied renormalization-group techniques to TMD PDFs, defined 
in a gauge-invariant way. We have shown by explicit calculation in the light-cone gauge 
of the one-loop gluon radiative corrections to the quantity fg/g {x, k±) that a contribution 
appears, which is proportional to Inp"*". This contribution gives rise to an anomalous dimen- 
sion that is formally equal to the universal cusp anomalous dimension and helps unravel a 
cusp obstruction in the composed gauge contour at light-cone infinity. The origin of this 
anomalous dimension can be traced to the renormalization effect on the junction point of the 
split contours, each associated with a gauge link and attached individually to a quark field — 
transforming it into a Mandelstam path-dependent field 49, 13]. Guided by this finding. 



we derived a generalized factorization rule for cusp- connected gauge links which contains an 
additional phase factor and worked it out. For gauge links joined along lightlike contours, 
this factor reduces to unity, while for more convoluted contours which run off to light-cone 
infinity in the transverse configuration space, this eikonal factor contributes an anomalous 
dimension due to the cusp. In this context, we emphasize that we verified that integrating 
over the transverse momenta in fg/g {x, k±) no artefact owing to the contour cusp remains, 
thus invigorating the validity of the standard integrated PDF. 

In order to eliminate the cusp anomalous dimension and recover the well-known results in 
a covariant gauge, where the gauge field vanishes at infinity, we proposed a new definition for 



the TMD PDF, which includes a soft counter term in the sense of Collins and Hautmann [43 



441 ]. in order to eliminate the contribution from the cusp anomalous dimension. This counter 
term enters in addition to the transverse gauge inks, previously introduced by Belitsky, Ji, 
and Yuan j33|, and comprises two eikonal factors caused by a particle-like current flowing 
along a cusped contour meandering from (0", — oo"^, 0_l) to (^^, oo+,^_|_) with a sharp bend 
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in the transverse direction. We have argued that each of these eikonal factors resembles 
in crucial aspects the "intrinsic Coulomb phase" found before by Jakob and Stefanis 



in a formulation of QED in terms of Mandelstam fields. In the present case, the cusp-like 
junction point of the two individual gauge contours plays a similar role as the so-called 
"particle behind the moon", postulated in QED in connection with the Lorentz-covariance 
restoration of its charged sector. Both quantities share the feature of being "hidden" at 
infinity and reveal themselves only in terms of (path-dependent) phases, being independent 
of external charge distributions (QED case) and unrelated to boundary conditions to avert 
light-cone singularities (TMD PDF case in QCD). The origin of the phenomenon is in both 
cases the same and peculiar to the inclusion of the individual path-dependent exponential 
into the field operator supposed to describe the quark as a Mandelstam field. No such 
effect appears in cases where the dynamics of the process allow one to use a direct contour 
between the two field points. In that case, one has to deal only with the connector which 
has well-known renormalization properties [3, 67 . 



The "intrinsic Coulomb phase" in QED tells us that each charged particle, though pri- 
mordially separated from its balancing counterpart, is still in harness with it. In the TMD 
PDF case, this phase accumulates effects due to the interaction of the struck quark with its 



target spectators, as pointed out by Ji and Yuan in [4l|] and reinforced by Belitsky, Ji, and 



Yuan in [33[ . However, the existence of a cusp at light-cone infinity went unnoticed, because 
in previous works the UV divergences of the TMD PDF were not considered. In |35| UV 
divergences were addressed, but only within the CoUins-Soper approach, which is formulated 
off-the-light-cone and, hence, the Inp"*" term does not appear there at all. 

The appearance of the cusp anomalous dimension in the present context is, in actual 
fact, not really surprising. We know from the so-called modified factorization of exclusive 
reactions that retaining transverse degrees of freedom amounts to the inclusion of Sudakov 
factors for each quark in the hard-scattering subprocess [7l[ — see for a review 

fli- The 

connection of the Sudakov factors to the cusp anomalous dimension within the modified 
factorization scheme was worked out in detail in [73] up to the level of the next-to-leading- 
order logarithmic accuracy. 

Our results may have a wide range of applications. Chief among them: 

• More precise data analyses of various experimental data on hard-scattering cross sec- 
tions. 

• Development of more accurate Monte-Carlo event generators (to estimate exclusive 



components of inclusive cross sections) [66|, |74 



Better description of polarized TMD PDFs and the phenomenology related to SSA 
and spin physics 0, i, M, 
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